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We include immunity against fire as a new parameter into 
the self-organized critical forest-fire model. When the im- 
munity assumes a critical value, clusters of burnt trees are 
identical to percolation clusters of random bond percolation. 
As long as the immunity is below its critical value, the asymp- 
totic critical exponents are those of the original self-organized 
critical model, i.e. the system performs a crossover from per- 
colation to self-organized criticality. We present a scaling the- 
ory and computer simulation results. 

PACS numbers: 05.40. +j, 05.70.Jk, 05.70.Ln 



Several extended noncquilibrium systems show scaling 
behavior over a wide range of parameter values and inde- 
pendently of the initial conditions. They are called self- 
organized critical (SOC) [Q and might explain the ubiq- 
uity of fractal structures in nature. By analogy to equi- 
librium critical phenomena, the question arises if SOC 
phenomena are universal, i.e. if the critical behavior 
depends only on few properties of the model as dimen- 
sion, conservation laws, or number of components. So 
far, examples for both universal and nonuniversal behav- 
ior have been found. While the critical behavior of the 
sandpile model Q seems to be robust with respect to 
various changes, the critical exponents of an earthquake 
model Q change continuously as function of a parameter 
which characterizes the degree of conservation. In equi- 
librium critical phenomena, such a continuous change of 
critical exponents occurs only in exceptional cases, e.g. 
the eight vertex model. Instead, critical behavior is char- 
acterized by a crossover between two fixed points, when 
certain parameters are changed. This crossover is de- 
scribed by scaling functions and a crossover exponent <fi. 
In this letter, we report for the first time a crossover phe- 
nomenon in SOC systems by including a new parameter, 
the immunity g, in the SOC forest-fire model ||. g is 
the probability that fire cannot pass from one tree to a 
neighboring tree. When the immunity assumes a criti- 
cal value g Cl the forest becomes very dense, and clusters 
that are burnt by a lightning stroke are identical to per- 
colation clusters. As long as the immunity is below its 
critical value, the asymptotic critical exponents of the 
model are those of the original SOC forest-fire model, 
i.e. the system performs a crossover from percolation to 
SOC. 

The forest-fire model is defined on a c?-dimensional hy- 
percubic lattice with L d sites. Each site is either occupied 
by a tree, a burning tree, or it is empty. The state of the 



system is parallely updated according to the following 
rules: 



1. Burning tree — > empty site. 

2. Tree — ► burning tree with probability 1 
n > 1 nearest neighbors are burning. 



if 



3. Tree — > burning tree with probability / if no near- 
est neighbor is burning. 

4. Empty site — ► tree with probability p. 

Starting with arbitrary initial conditions, after some time 
the system approaches a steady state the properties of 
which depend on the parameter values but not on the 
initial state. We always chose the lattice size so large 
that the steady state is also independent of the boundary 
conditions. Let p t be the density of trees, p e the density 
of empty sites, and pj the density of burning trees in the 
steady state. Large forests and therefore large fires occur 
only if the lightning probability / is much smaller than 
the tree growth probability p. If additionally large and 
small fires look alike, the system becomes scale invari- 
ant. Small fires live only for few time steps and are ex- 
tinguished before new trees grow in their neighborhood. 
Scale invariance therefore can only be observed if tree 
growth is so slow that large fires also are extinguished 
before new trees grow at the edge of the burning forests. 
We conclude that SOC behavior occurs in the forest-fire 
model if the following conditions are satisfied: Lightning 
occurs seldom compared to tree growth, and tree growth 
is much slower than the lifetime of a fire, i.e. 



/ < p < (f/pf 



(1) 



with an appropriate exponent v' || . Eq. (Q) represents a 
double separation of timescales. Since in the steady state 
the mean number of burning trees equals the mean num- 
ber of growing trees, the mean number of trees destroyed 
by a lightning stroke is 
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In the last step, we have neglected the fire density pf 
which vanishes in the limit of perfect time scale sepa- 
ration. In d > 2 dimensions and for values of the im- 
munity below its critical value, the critical forest density 
pi = lim j /p^o Pt is smaller than 1, and Eq. (^|) repre- 
sents a power law s oc (//p)~ 7 with 7 = 1, indicating a 
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critical point at f/p = 0. Close to this critical point, i.e. 
if / <C p, there is scaling over many orders of magnitude. 

Several critical exponents characterizing this scaling 
behavior have been defined || 1| . The forest density sat- 
isfies for small //pa power law 



Pt 



Pt oc (f/ P ) 1/S , 



(3) 



as first stated in g and |6J]. 

Let s be the number of trees destroyed by a fire. The 
size distribution of fires is a power law fel 



with 



sn(s) oc s 1 T C(s/s max ) 



« (f/pY 



(4) 



(5) 



Inserting Eq. (g) and Eq. (||) into Eq. (^), we obtain 
the scaling relation A = 1/(3 — r), which is valid as long 
as p% < 1. In general, if we allow also the case pi = 1, as 
relevant for g = g c , s is given by a power law s oc (//p)~ 7 , 
with a value for 7 which may be different from 1 , and we 
have the more general scaling relation 



A = 7 /(3-r) 



(6) 



We also define the radius R(s) of a forest cluster that 
has just been burnt down. Its fractal dimension p is given 
by 



(7) 



(8) 



R(s) oc s 1 ^. 
The correlation length £ is defined by 

Together with Eqs. (§), (§) and (0), this gives 

£oc(//p)- I/ with l / = A/ M . (9) 

In the SOC forest-fire model without immunity, these 
exponents have already been determined by computer 
simulations In d — 2 dimensions, they are f| 

r soc = 2.14(3), A soc = 1.15(3), p soc = 1.96(1), 
i/ soc = 0.58, l/S soc = 0.48(2). (10) 

In d — 1 dimension, their values have been derived an- 
alytically M. In ||, it is also shown that the critical 
exponents in two dimensions do not change when the 
lattice symmetry is changed or next-nearest neighbor in- 
teraction is included. 

We now consider the SOC forest-fire model with non- 
vanishing immunity, i.e. for g > 0. Immunity was first 
introduced in Q into a version of the forest-fire model 
g§ which is not SOC This model shows a 

percolation-like phase transition when the immunity ap- 
proaches a critical value g c . In the present paper, the 



immunity is for the first time included into the SOC 
forest-fire model. It is defined differently than before: 
While in || the immunity is a property of trees, in the 
present paper it is a property of bonds between neighbor- 
ing trees (see rule 2. above). This has two advantages: 
The simulation program is less complicated, and the crit- 
ical immunity is known to be exactly g c — 0.5 which is 
just 1 minus the percolation threshold for bond percola- 
tion. 

When the immunity is different from 0, not all trees 
that are neighbors of a burning tree catch fire, and con- 
sequently the fire does no longer burn forest clusters 
but clusters of trees that are connected by non immune 
bonds. With increasing immunity, the forest density in- 
creases, since fewer trees are burnt. At the critical immu- 
nity <7 C , the critical forest density is pi = 1. Then we have 
the following situation: The forest is completely dense in 
the limit f/p — > 0, and clusters that are destroyed by 
fire are percolation clusters of bond percolation. Conse- 
quently the exponents r and p are given by percolation 
theory 



t(9c) 



' pcrc 



187/91 ~ 2.05 



and 



p(g c ) = p c = p peIC = 91/48 ~ 1.90. 



(11) 



(12) 



Bonds that are immune during one time step might be 
non immune during the next one, and consequently light- 
ning strokes at the same site burn down different clusters 
at different times. When f/p is finite, there is a cutoff in 
cluster size, since large fires are stopped by empty sites 
that have been left from earlier fires. The mean forest 
density is no longer 1. We determined the critical expo- 
nents A, 5, and v at g — g c by computer simulations in 
d = 2 dimensions and obtained 

A c = 0.92(3), l/6 c = 0.15(1), v c = 0.484(2). (13) 

The scaling relations Eq. (JsJ) and Eq. @ still hold at 
g = g c and are confirmed by our simulations. Using Eqs. 
(0) and (0), we obtain 



S oc (f/p)~^ with 7 c = 1 - l/Sc 



(14) 



Our simulations yield "f c = 0.84(2), in agreement with 
Eq. @. 

When the immunity is just below its critical value 
((.9c — g) *C 1), the situation becomes more complicated. 
On small length scales, a system close to the percolation 
threshold cannot be distinguished from a system exactly 
at the percolation threshold. On large length scales, how- 
ever, the difference can be seen. If the initial state were a 
completely dense forest, a finite fraction of all trees were 
connected by non immune bonds for g < g c . This in- 
finite cluster would soon be destroyed by lightning and 
would never occur again due to Eq. (||). Consequently 
the critical forest density is pi < 1 for g < g c . Imag- 
ine for a moment that bonds are permanently immune 
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or non immune (in reality, rule 2. says that they might 
be immune during one time step and non immune dur- 
ing the next one) . Then the dynamics on the sites which 
initially belonged to the infinite cluster are completely 
decoupled from dynamics on the remaining sites. We can 
consider the sites which belonged to the infinite cluster 
as an independent subsystem where no immunity exists. 
They form a two-dimensional lattice with another sym- 
metry than the original hypercubic one. Since we already 
know from earlier simulations that the exponents of the 
SOC forest-fire model do not change when lattice sym- 
metry is changed, we expect that the critical exponents 
on our subsystem are just the ones which we obtained 
for g = and which have been assigned an index " soc" . 
The remaining sites are connected by non immune bonds 
to form finite clusters of bond percolation. On each of 
these clusters, the dynamics are independent of the other 
clusters and of the infinite subsystem. The largest of 
these clusters have a radius of the order of the percola- 
tion correlation length £ pcrc oc (g c — g)~ u v mc . When f/p 
becomes very small, these finite clusters are very rarely 
struck by lightning, and consequently the mean tree den- 
sity on these clusters approaches the value 1 in the limit 
f/p — > 0. Each time lightning strikes such a cluster, 
there is a fire which has the size and the fractal dimen- 
sion of a finite percolation cluster. As long as f/p is so 
small that the correlation length £ oc (/ /p)~ v " is much 
larger than the percolation correlation length £ p0 rc, ah 
large fires occur on the infinite subsystem, and the expo- 
nents v, 5, and A are those of the subsystem, i.e. the SOC 
ones. The exponents r and \i of the large fires are by the 
same reasoning also the SOC ones. But for fires with a 
radius smaller than the percolation correlation length, r 
and p assume their percolation values since the infinite 
subsystem contains only a small portion of all sites and 
consequently most of the small fires occur on the finite 
clusters. When f/p becomes so large that the forest clus- 
ters on our infinite subsystem are no more larger than the 
percolation correlation length, our system is dominated 
by the dynamics on the finite clusters and becomes in- 
distinguishable from a system at g = g c . The exponents 
v, S, and A then are identical to those at g = g c , and 
the exponents r and ji are those of percolation theory on 
all length scales up to the correlation length. So far, we 
considered the case of small g c — g. When the immunity 
is far away from its critical value, the infinite subsystem 
contains a large portion of all sites, and the percolation- 
dominated behavior cannot occur any more. 

Unfortunately, all these considerations are based on 
the assumption that bonds are permanently immune or 
non immune which in reality is not the case. Conse- 
quently, there exists no subsystem which is decoupled 
from the rest of the system. Nevertheless, the main con- 
clusions should remain valid: On length scales smaller 
than the percolation correlation length, the system can- 
not be distinguished from a system at g — g c . When f/p 
becomes very small, there are fires which spread further 
than the percolation correlation length. These fires are 



stopped by empty sites that were created by earlier fires. 
This is again the same mechanism as in the limit g = 
or on the infinite subsystem in the case of fixed immune 
bonds: fires that would spread indefinitely if there were 
no empty sites are stopped by empty sites. We conclude 
that these large fires lead again to the critical exponents 
Asoc; ^soci and 5 soc . 

Led by these considerations, we make the following 
scaling ansatz for the correlation length: 

It is plausible that the crossover from percolation-like to 
SOC behavior takes place when f/p becomes so small 
that the correlation length exceeds the percolation corre- 
lation length, which suggests that the crossover exponent 
<f> is 

4> = v c /v perc . (16) 

The scaling function F(x) is constant for small x and is 
fx x {-fsoo-v c }/4> f or i ar g e x Analogous scaling laws hold 
for s max and p c t — p t . We already mentioned above that 
the critical forest density is p c t = 1 at g c . We therefore 
expect an additional power law 

l-p c t (g)K(g c -g)y. (17) 

The exponent y is obtained from the scaling ansatz 

W t = (//P)^(§/^). (18) 

In the limit f/p — » 0, the forest density becomes inde- 
pendent of f/p and assumes a value p% ^ 1. Therefore 
G(x) oc x x l^ c for large x, yielding 

V = Vperc/v c 5c- (19) 

Our simulations confirm all these results. They were 
performed using the same method as in |,§ for lat- 
tices of up to 8192 2 sites and values of f/p down to 
1CP 6 . The immunity varied between (g c — g) = 0.02 and 
(g c — g) — 0.001. The range of f/p is limited by finite-size 
effects for small f/p and noncritical behavior for large 
f/p, the range of g was selected from the condition that 
at least part of the crossover region is covered. Fig. |l| 
shows the scaling function for the correlation length F(x) 
for different values of g c — g. The scaling ansatz Eq. ( jig ) 
is well confirmed since all curves coincide. The dashed 
line represents F(0) as obtained from the simulations at 
g c . We also checked the scaling relation Eq. (|l9|). Fig. || 
shows the critical forest density as function of g c . We ob- 
tain y = 0.43(2), in agreement with Eq. ( |l9| ) (remember 
that Vpcrc = 4/3). Since the difference between t soc and 
Tpcrci as well as between psoc and Pperc 

is very small, the 

crossover behavior of n(s) and R(s) could not be evalu- 
ated. 
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In §, we also denned exponents u 1 , fi' and a describ- . FIG - 2 - The critical forest densit y as Unction of the immu- 
ing the temporal behavior of the fires. The crossover in nlty ' 
// and v' , which enters the condition for time scale sep- 
aration Eq. (^), is analogous to the crossover in fj, and v 
. The change in a is too small for any crossover to be 
observable. 

To conclude, we have shown by analytic arguments 
and by computer simulations that the forest-fire model 
performs a crossover from percolation to SOC when the 
immunity is close to its critical value. This crossover is 
characterized by scaling functions which are defined in 
the same way as in crossover phenomena at equilibrium 
phase transitions. 

Although all simulations were performed in d = 2 di- 
mensions, we expect that this crossover behavior can also 
be observed in higher dimensions. In d = 1, the critical 
immunity is g c = 0, and no crossover can take place. For 
d > 6, simulations suggest that the critical exponents as- 
sume their mean-field values which are identical to those 
of percolation Consequently there is no crossover 

in d > 6 dimensions. 
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FIG. 1. Crossover scaling function F(x) for the correlation 
length for different values of the immunity. The dashed line 
represents F(0) as obtained at g = g c . 
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